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Abstract 

Two classical results of Stafford say that every (left) ideal of the n-th 
Weyl algebra A n can be generated by two elements, and every holonomic 
^-module is cyclic, i.e. generated by one element. We modify Stafford's 
original proofs to make the algorithmic computation of these generators 
possible. 

1 Introduction 

Let k is a field of characteristic 0, and A n = A n (k) = k{x\, ...,x n ,di, ...,d n ) be 
the n-th Weyl algebra, which is an associative fc-algebra generated by x's and 
d's with the relations diXi — Xidi + 1 for all i. This algebra may be thought of 
as the algebra of linear differential operators with polynomial coefficients. 

There are several things that are nice about the Weyl algebra. First of all the 
dimension theory can be developed for it; this is done, for example, in Chapter 
1 of Bjork It is shown that the Gelfand-Kirillov dimension of A n equals 2n, 
moreover, if M is a nontrivial A n -module, then n < dimM < In. The modules 
of dimension n (minimal possible dimension) constitute the Bernstein class. 

One of the distinctive properties of the modules in Bernstein class, which 
are also called holonomic modules, is their finite length. Below we shall show 
that this property implies that every holonomic module can be generated 
by one element. 

Another striking fact, which is very simple to state, but quite hard to prove, 
is that for every left ideal of A n there exist 2 elements that generate 
it. 

Both statements were proved by Stafford in M, also these results appear in 
. Unfortunately, the arguments given by Stafford can't be converted to algo- 
rithms straightforwardly. There are several obstacles to this, many of which one 
can overcome with the theory of Grobner bases for Weyl algebras. However, the 
main difficulty is that both proofs contain an operation of taking an irreducible 
submodule of an A n -module. To our best knowledge, there doesn't exist an al- 
gorithm for this; moreover, even if such algorithm is invented one should expect 
it to be quite involved. 



We were able to modify the original proofs in such a way that computations 
are possible and implemented the corresponding algorithms in the computer 
algebra system Macaulay 2 B . 

We have to mention that in their recent paper || Hillebrand and Schmale 
construct another effective modification of Stafford's proof which leads to an 
algorithm. We shall discuss the differences of their and our approaches in the 
last section. 

2 Notation Table 

For the convenience of the reader we provide the notation lookup table. All of 
the symbols listed below show up sooner or later in the paper along with more 
detailed definitions. 



k 


is 


a (commutative) field of characteristic 0, 






A r (k) = k(xi, ...,x r ,d%, ...,d r ), 


A 


is 


a simple ring of infinite length as a left module over itself, 


D 


is 


a skew field of characteristic 0, 


K 


is 


a commutative subfield of D, 


S 




D(x)(d), 


S (m) 




Sei + ... + Se m , a free S'-module of rank m, 


i j • • • j 8 m 


is 


a finite set of X-linearly independent elements in K(x,d), 






rn 

J2 aS if £ i e ^ (m) , (aeSJe K(x,d)), 

i=l 






Sa(aJ), ideal of 5 (m) , 


V r 


is 


the quotient ring of A r , 


Kr 




V r (x r+1 , ...,x n )(d r+ x, ...,d n ), 


*-V 




T> r (x r+ i, ...,x n )(d r+ i). 



With exception of some minor changes we tried to stick to the notation in 



3 Preliminaries 

Several useful properties of Weyl algebras are discussed in this section. Also, 
we introduce a few rings that will come handy later on. 

3.1 A n is simple 

To see that A n is simple, i.e. has no nontrivial two-sided ideals, we notice that, 
for / = J2i x ai d l3i E A n \ {0} in the standard form, df jdx T = d r f — fd r for 
r = l,...,n, where df/dx r is the formal derivative of the above expression of 
/ with respect to x r . Similarly, df /dd r — fx r — x r f for the formal derivative 
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with respect to d r . Note that these formal derivatives as well as all the multiple 
derivatives of / belong to the two-sided ideal A n fA n . 

Now assume x a d /3 is the leading term of / with respect to some total degree 
monomial ordering. We are going to perform \a\ + \{3\ differentiations: for all 
i = 1, ...,n differentiate / oq times with respect to x, and /Jj times with respect 
to d%. Under such operation the leading term becomes equal to Yli=i a i W an( A 
all the other terms vanish. Since the derivatives of / don't leave A n fA n , we 
showed that there is a simple algorithm to find such Sj , r, G A n that 

m 

^sjr, = 1. 

i=l 

Hence, A n fA n = A n , so A n is simple. 

3.2 An is an Ore domain 

Proposition 1 A n is an Ore domain, i.e. A n f D A n g ^ and fA n n gA n ^ 

for every f,g£ A n \ {0}. 

Proof. See the proof of Proposition 8.4 in Bjork H. ■ 

Let us point out that using Grobner bases methods (see next subsection) we 
can find a left (right) common multiple of /, g £ A n \ {0}, in other words we can 
find a nontrivial solution to the equations af = bg and fa = gb where a and b 
are unknowns. 

3.3 Grobner bases in A n 

As we mentioned before, the notion of Grobner basis of a (left) ideal can be de- 
fined for Weyl algebras in the same way as it is defined in the case of polynomials. 
Moreover, Buchberger algorithm for computing Grobner bases works, leading to 
algorithms for computing intersections of ideals, kernels of maps, syzygy mod- 
ules, etc. A good reference on Grobner bases for algebras of solvable type is 

3.4 More rings 

There is a quotient ring D associated to every Ore domain A. Ring D is a skew 
field that can be constructed both as the ring of left fractions a~ 1 b and as the 
ring of right fractions cd _1 , where a, 6, c, d G A. There is a detailed treatment 
of this issue in [|J . 

Let D be a skew field, we will be interested in the ring 5* = D(x){d), which 
is a ring of differential operators with coefficients in D{x). It is easy to see that 
S is simple. 

Since the Weyl algebra A r is an Ore domain, we can form its quotient 
ring, which we denote by T> r . The S we are going to play with is S r = 
T> r (x r+ 2, x n )(x r+ i)(<9 r+ i). Let us state without proof a proposition which 
shall help us to compute Grobner bases in S r . 
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Proposition 2 Let F — {/i, ...,fk} C A„ is a generating set of left ideal I of 
S r . Compute a Grobner basis G = {gi, g m } of A n ■ F with respect to any 
monomial ordering eliminating d r+ \. Then G is contained in S r (~l A n and is a 
Grobner basis of I. 

4 Holonomic modules are cyclic 

In this section we consider a simple ring A such that A has finite length as a 
left module over itself. Note that A n is such a ring. 

Theorem 3 Every left A-module M of finite length is cyclic. In particular 
every holonomic A n module is cyclic. 

Suppose we know how to compute a cyclic generator for every module M' 
of length less than I. For length such generator would be 0. 

Consider a module M of length I. Take / a € M. If M = Act then we are 
done. If not then since l(M/Aa) < I by induction we can find (3 such that its 
image in M/Aa is a cyclic generator. Now M = A ■ {a, (3} and what we need 
to prove is 

Lemma 4 Let M be a left A-module of finite length and a, (3 S M. Then there 
exists 7 G M such that Aj = Aa + Af3. 

Proof. Define two functions l\ and l 2 for pair (a, (3). 

h(a,0) = length^/?) 

h (a, [3) = length ( {Aa + A(3) /Aa) . 

Let also introduce an order < on the set of pairs (a, (3) e M x M: 

{a',[3')<{a,(3) ^ (h(a\ /?'), h(a', (3 1 )) < lex (h(a, (3), l 2 (a, (3)) 
«■ h(a' , (3') < h(a, f3) 

OR {h(a',(3') = h(a,(3) AND h{a',0) < h(a,0)) 

Suppose for any pair (a', [3') < (a, (3), we can find 7' 6 M such that Aj' — 
A-{a',f3'}. 

Let the ideals L(a) and L{(3) in A be the annihilators of a and (3 respectively. 
Since lcngth(^4) = 00, we know that L(a) ^ 0; pick any element 0^/6 L(a). 
Since A is simple we can find Sj, G A, I = 1, M such that 

m 

J>i/ri = l. (1) 

i=l 

Consider two cases: 

1. There is some r — r t such that L([3) + L(a)r = A. 
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2. The opposite is true. 

Case 1. We can write 1 = E a r + Eg for some E a , Ep 6 A such that E a a = 
and Epfi — 0. Let 7 = a + r/3. 
Now we can get /3 from 7: 

- (£ Q r + Ep) /3 = E a r(3 = E a a + E a rj3 = E al . 

Hence (3 6 A7 and since a — 7 — r/3 the module M = Aa + A(3 is indeed 
generated by 7. 

Case 2. From @j it follows that £ L(/3) + A/r, = A, hence, £ A(/r 4 /3) = 
Aft, so there is r = r% such that 

A(/r/3) g Aa. (2) 

Since we are not in case 1, L((3) + Afr C + L(a)r ^ A. Take this 

modulo L(/3) to get 

A(frft) S + A/r) /L(0) C A/L(0) S A/3, (3) 

so A(frf3) is proper in A/3. 

The last statement implies li(a,fr(3) < l\{a,(5), hence, (a,fr(3) < (a, (3), 
so by induction hypothesis we can find 7' 6 M such that A7' = A(/r/3) + Aa. 

Now @ guarantees that foil', P) < h{&, /3), and by induction we can find 7 
for which 

A7 = Ay' + A/3 = A(/r/3) + Aa + A/3 = Aa + A/3. 



Remark 5 There is an algorithm that finds a cyclic generator for a holonomic 
left module over a Weyl algebra, since every step in the proof of the Lemma Q 
is computable. The most non-trivial and time consuming operation is producing 
the annihilators L(a + r/3) and L(fr/3) in the proof of Lemma [J provided L(a) 
and L{j3). This is done using Grobner bases technique. 

We have programmed the algorithm corresponding to the proof of Theorem 
[| using Macaulay 2. 

Example. Let us view the ring of polynomials k[x] as an Ai-module un- 
der the natural action of differential operators. It has an irreducible module, 
because starting with a nonzero polynomial / we can obtain a nonzero con- 
stant by differentiating it deg(/) times. The module M = k[x] 3 is the direct 
sum of 3 copies of k[x], is holonomic (length(M) = 3) and is generated by 
vectors (1, 0, 0), (0, 1, 0), (0, 0, 1). Our algorithm produces a cyclic generator 
7 = (x 2 , x,l) and its Ai-annihilator £(7) = Aid 3 . 
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5 Ideals are 2-generated 



In this section we give an effective proof of 

Theorem 6 Every left ideal of the Weyl algebra A n can be generated by two 
elements. 

Proof for A\. In this case the theorem follows from the fact that module 
Ai/J is holonomic for any nonzero ideal J of A\. 

Indeed, let I be a left ideal of A\. Pick / G I and set J = Aif . Then I/J is 
a submodule of the holonomic module A\j J, hence, is holonomic. By Theorem 
H there is g G I/J such that A\g = I/J. Find a lifting g G Ai such that g = g 
mod J. Elements / and g generate /. ■ 

However, the theorem for n > 1 makes a much tougher challenge. 
5.1 Lemmas for S 

Let us explore some properties of S = D(x)(d), the ring of linear differential 
operators with coefficients in rational expressions in x over a skew field D. 

Let K be a commutative subfield of D, let 6i,...,S m be a finite set of K- 
linearly independent elements in K(x,d) C S, and let S( m ' = Sei + ... + Se m 
be a free 5-module of rank m. 

Also define cr(a, f) G to be the following sum a(a, f) — YliLi a ^if £ i^ 
and P(a, f) — Sa(a, f) the submodule of S ,( - m - > generated by er(a, /). Note that 
o~(a, f) is S'-linear in a and respects addition in /. 

Lemma 7 Let ^ a G S and let M be an S-submodule of S^ m ^ generated by 
W(aJ)\f G K{x,d)}. Then M = S^. 

Proof. Without loss of generality let us assume that a G D(x,d): if not wc 
can always find such p G D[x] that pa G D(x,d). 

Fix a monomial ordering that respects the total degree in x and d. For 
vector v = J2 v i £ i e {D{x,d))^ denote by lm(?;) the largest of the the leading 
monomials of the components i>; of v in this ordering. 

Now start with vector v = = o~(a, 1); its components Vi = aSi are D- 
linearly independent. Note that computing expressions n(v) — dv — vd and 
x(v) = vx — xv has an effect of differentiating each component of v formally 
with respect to x and d respectively. These operations lower the total degree of 
v by 1 if the differentiation is done with respect to a variable that is present in 
\m(v). Also, it is not hard to see that they keep us in module M; for example, 
for vq we have tt(vq) — Ovq — vq8 = da{a, 1) — a(a, d). 

Run the following algorithm: initialize v := vq, while lm(w) contains an x 
set v := 7r(f), then while \m(v) contains a d we set v :— x( v )- Since each step 
lowers the total degree of v by 1, this procedure terminates producing vector 
w G M of total degree 0. 
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Hence, w = Wi 1 Ei 1 + ... + Wi t Si t where ^ Wi j € D for j = 1, ...,t. Via 
multiplying on the left by the inverse of we can get the relation 

e il = a 2 e l2 + ... + a t e it (4) 

with ctj G D for j = 2, t. 

Now take and reduce it using . We get vector whose ii-th com- 
ponent is and the remaining components are £)-linearly independent, since 
the components of v^ ' are. 

Repeat the above algorithm for v = and so on. At the end we get a 
vector which is a scalar multiple of Si for some i, hence e, £ M. Using relations 
we see that all basis vectors Ej, for j = 1, m, are in M. m 

Remark 8 From the proof it follows that given a submodule M of S^ m ' and 
a G S one can find f £ K(x, d) such that o~(a, f) £ M algorithmically. 

The next lemma is central in the proof of the result. Note that every step 
of the proof of the lemma can be carried out algorithmically. 

Lemma 9 Let M be an S-submodule of S^ m ' — Sei + ... + Se m such that 
length(S^ /M) < oo. We can find f e K{x, d) such that S (m *> = M + P{a, /). 

Proof. Let / = length(5 1 ' m > /M). Assume the assertion is proved for all M' 
such that length(5( m V M ') < l - Remark | says that we can find an / 6 K(x, d) 
such that a (a, f) doesn't belong to M. 

For t G S, g S K(x, d) let us define two S'-modules 

Ni = M + Pi, where Pi = P(a, g), 
N 2 = M + P 2 , where P 2 = P(ta, g). 

Claim. There is a module M' such that M C M' C M + P(a,f), t G S, 
and g G K(x,d) for which 

ta(a, f) G M, 

M' + P(ta, g) = S (m \ 

Ni = N 2 . 

To prove this we employ (second) induction on length(M'/M). We start 
with M' = M + P(a,f). We can find such that taJ^Stfs, e M ; 

it follows from S being Ore. By the first induction hypothesis, for M' and ta 
there exists g G K{x, 0} such that W + P(ta, g) = 5 (m) . Notice that Ni D N 2 
and M' + P t = 5 (m) for i = 1,2. Also for i = 1, 2 we have 

SW/Ni = {M 1 + Pi)/(M + Pi) = M'/(M + M' n P l ). 

If length(5 (m) /^i) = lcngth(S' (m) /A r 2) then N x = N 2 and we are done. We 
are done as well if N\ — S^ m \ If both conditions above fail, by looking at the 
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right hand side of 5.1 we determine that M" = M + M' n P\ both contains 
M and is contained in M' properly, plus length(M"/M) < length(M'/M). Set 
M' := M" and repeat the above procedure. 

To finish the proof of the lemma we take M ', t, g as in the claim and assert 
that N' — M + P(a,f + g) equals S^. Indeed, a(taj + g) = ta(a, f) + 
a(ta,g) = a(ta,g) modulo M, so N 2 C N'. But Ni = N 2 , thus a(a,g) € N', 
hence, a(a,f) = a (a, f + g) — a(a ) g) 6 N' . Now we see that M' C N' and 
P 2 C N'. Since M' + P 2 = S (m) , we proved N' = S^ m \ m 



5.2 Lemmas for lZ r 

At this stage we shall specify the components in the definition of S = D(x)(d). 
We set D = V r (x r+2 , x n ), x = x r+ i and d = 9 r +i, so that new S is equal to 
<S r = T> r (x r +i, x r + 2 , x n ){d r +\) which is a subring of lZ r . Also the commuta- 
tive subfield K of D that showed up before is replaced by the k, the coefficient 
field from the definition of A n = A n (k). 

Proposition 10 Let 5\, S m be a finite set of K -linearly independent elements 
in K(x r+ i, d r+ i) and let ^ p 6 A r+ i[x r+2 , x n ]. Let 5 , ' m+1 ' = Seq + Se± + 
... + Se m be a free S-module of rank m + 1 And let S'( m + 1 ) / 9 C S^ m+1) be its S- 
submodule generated by {pSQ,pe\, ...,pe 2 }. Then there exists some f £ K such 
that 

S {m+1) = S (m+l) p + gfa + ^ + + 5m f £m) . 

Proof. Follows from Lemma || ■ 

Lemma 11 Let q £ A r [x r +i, x n ] and let a\, at be a finite set in A n . 

Then there exists some ^ p £ A r [x,+i, ...,x n ] such that paj £ A n q for all 

3- 

Proof. See the proof of Lemma 8.5 in Bjork H. ■ 

Let us point out that once we know the statement of the lemma is true, we 
can compute the required p by finding a Grobner basis of the module of syzygies 
of the columns of the matrix 

/ at q ... \ 
a 2 q ... 

\ a f ... q ) 

with respect to a monomial order that eliminates d r +i, d n and such that 
£\ > e 2 > ... > Et+i where £1,62, ■■■,£t+i is the basis (£j corresponds to the z-th 
column) of the free module A^ 1 containing our submodule of syzygies. Such 
Grobner basis is guaranteed (by Lemma to contain some syzygy producing 
the relation pe\ + b 2 s 2 + ... + b t+ iet+\ = where p £ A r [x r +i, x n ], bi £ A n 
for i = 2, n. It is not hard to see that this is the p we need. 
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Lemma 12 Let =/= q € A r+ i[.x r+ 2, ■ a^d 1st u,v 6 urai/i a / 0. TTien 

£/iere «s some / £ A„ sucft £/iai 7?. r = 7?. r g + lZ r (u + u/). 

Proof. Consider the following subring of A n obtained by "removing" x r +\ and 

dr+l- 

Af^- = k(xi, ...x r ,x r+2 , ...,x n ,di, ...,d r ,d r + 2, ...,d n ). 

Now A n — A^p^ <£>k k(x r+ i,d r +i), so we can write v — 6\gi + ... + S m g m 
where 5i,...,S m are elements of k(x r +i, d r +i) linearly independent over k and 
gi, ...,g m G ^4^rjrp The ring A^^ is simple, since it is a Weyl algebra, thus we 
can find such hi,..., hi G ^rpj that 

m / 
i=l j=0 

Since ^rpj is a subring of 7^.,- it means that lZ r = YlYl !Z r gihj. 

Sublemma. For any bi, b m G ^rpj there exists some f G fc{av+i, <9 r +i) 
such that 

TZ r q + TZ r u + TZ r bi + ... + 1Z r b m = TZ r q + TZ r (u + difbi + ... + 6 m fb m ). 

Proof. It follows from Lemma [TI] that there is =/= p G A r [x r+ i, ...,x n ] such 
that pbi, —,pb m G A n q as well as pu G A n q. With the help from Proposition 
[tO]wc get / e k{x r+ i,d r+ i) such that S ( - m+r > = S {m+ ^ p + S(e Q + 5 1 fe 1 + ... + 
SmfEm) and since S is a subring of lZ r we have 

T#»+ 1 ) = 7^ m+1 V + K r (e + 6if £l + ... + S m fe m ). (5) 

Now map £o i— » u and Si \— * bi for all i\ this map from 7^™ to lZ r has its image 
equal to 1Z r q + lZ r u + 1Z r bi + ... + lZ r b m and maps the right hand side of (||) to 
a subset of lZ r q + 1Z r (u + difbi + ... + S m fb m ), because pu, pbi, pb m G A n q. 
Moreover these two expressions are equal, since it is easy to see that the latter 
is contained in the former as well. ■ 

Proof of lemma continued. We apply our Sublemma to bi = g t h\ (i = 
1, ...,to) to get fi G k(x r +i, d r +\) such that 

m rn 
lZ r q + 7Z r u + lZ r gihi = lZ r q + 7Z r (u + Sifigihi). 
3 = 1 3=1 

Since v = Sigi + ... + S m g m and since fi commutes with all gi, the last equation 
transforms into 

lZ r q + lZ r u + lZ r gihi — lZ r q + !Z r (u + vfihi). 

Now reapply the Sublemma with u replaced by u + vfihi and bi = gihi 
(i = 1, to). As in the first step we get 

lZ r q + lZ r u + lZ r gihi + 'R- r gih2 

= TZ r q + TZ r (u + Vflhl) + ^2 ^-rffi^2 

= lZ r q + lZ r (u + vfihi + fih-i) 
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for some f% S fc(x r +i, d r +i)- After / many steps we arrive at 

7Z r = lZ r q + 7Z r u + 7Z r gihj = lZ r q + lZ r (u + v fihi), 

»=1 J=l 3=1 

which proves the lemma with / = + ... + fihi. ■ 
The following lemma follows from the previous one. 

Lemma 13 Let < r < n — 1 and let ^ q G v4 r +i[:r r +2, ain] fei 
u,v <E A n with v =/= 0. TTien there is some f s A n ,q' G ^4 r [a; r +i, ...,£„] swc/i t/ia£ 
g' e A„ g + A n (u + ?;/). 

Proof. It is easy to see that this lemma is equivalent to the previous one. ■ 
5.3 Final chords 

Proposition 14 (r) Let < r < n, there is some q r G A r [x r+ \, x n ] and 
d r , e r G A n such that q r c G A n (a + d r c) + A n (b + e r c). 

Proof. The statement is true for r = n, since A n is Ore and A n cC\{A n a + A n b). 

Fix r. Assume that the statement is true for r + l,...,n, then there exist 
q r +i, dr+i, e r +i such that q r +ic G A n a' + A n b', where a' = a + d r+ \c and 
b' = b + e-r+ic Hence we can write q r +\c — h\a! + h 2 b' , where we can take 
h\hi ^ since A n a' D A n b' ^ 0. Also since h\A n n h2A n ^ we can also find 
gi, g 2 satisfying hig% + h 2 g 2 = 0, and since A n q r+ icnA n b' ^ there are s, t such 
that sq r +ic = tb' . Using Lemma |l3| to q = q r +i with u = and v = tg2, we 
get q r — q' and / such that q r = p\q r +\ + p 2 tg 2 f for some pi,p 2 - Summarizing, 
there exist such hi, h 2 , gi, 32, s, t,pi,p 2 G A n \ {0} that 

q r =piq r +i +P2tg 2 f 
q r +ic = hia' + h 2 b' 
higi + h 2 g 2 = 
sq r+ ic = tb' 

Using these 4 equations, make the following calculation: (In each section the 
underlined terms sum up to 0.) 

q r c = piq r +ic + p 2 tg 2 fc 

= piq r+ ic - p 2 sq r+ ic 

+ P2tg 2 fc + p 2 tb' 

= (Pi ~ P2s)q r +ic + p 2 t(b' + g 2 fc) 

= (Pi ~ P2s){hia' + h 2 b') + p 2 t(b' + g 2 fc) 

= (Pi -p 2 s)hia' + {pi -p 2 s)higifc 

+ (Pi - P2s)h 2 b' + (pi - p 2 s)h 2 g 2 fc + p 2 t(b' + g 2 fc) 

= (Pi - P2s)hi(a' + gifc) + {{pi - p 2 s)h 2 + p 2 t)(b' + g 2 fc). 
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Thus, with d r = d r+ i + gifc and e r = e r +\ + g 2 fc the conclusion of the 
proposition holds. ■ 

The proposition above (for r = 0) shows that by "elimination" of variables 
di one at a time we can get such d, e S A„ that q$c G A„(a + cfc) + A n (fr + ec) 
where go €E fc[xi, x n \. This proves a 50% version of Theorem y: 

Theorem 15 Every ideal of k{x\, x n ){di, d n ) can be generated by two 
elements. 

To go other 50% of the way one has to do a similar kind of "elimination" 
of Xi-s. This amounts to making copies of all lemmas that we stated for a 
slightly different set of rings. The trickiest part is considering ring S r ' = 
k(x%, x r )(av+i, d r +\) instead of S r . In other words instead of a ring of type 
D{x)(d) where D is a skew field, we have to consider the first Weyl algebra 
Ax(K) where K, is a (commutative) field. Fortunately, analogues of Lemmas [?] 
and y for the latter ring can be effectively proved along the same lines. 
Examples. (1) Consider A3. For a = di,b — &2 : c = 83 one can show that 
A3 ■ {a, b, c} = A3 • {a, b + x\c\. Indeed, the following calculation displays it: 

c = (— x\dz - d 2 )a + di(b + xic). 

(2) Another example is produced by our algorithm implemented in Macaulay 
2. Let a = di + X3, b = 9| + x% + x|, c = 83 + x\. Then the ideal A 3 • {a, b, c} 
is generated by d\ + 23 and d\ + (x\xz + xi)d 3 + x^x^ + x\ + x\ + x 2 - 

(3) In case of A\ we can construct a more efficient algorithm based on the 
proof of Theorem ^| given for this special case. Here is a Macaulay 2 script 
computing 2 generators for the annihilating ideal / C Ai(Q) of the set of poly- 
nomials {ax 4 + bx 6 + cx 8 + dx 10 a, b, c, d e Q} C Q[x}. 

il : load "D-modules .m2" ; load "staf f ord .m2" ; 

13 : R = QQ[x,D, WeylAlgebra=>{x=>D>] ; 

14 : L = {4,6,8,10}; 

15 : I = ideal gens gb intersect apply(L, i->PolyAnn x~i) ; 

4 4 3 3 2 2 11 

o5 = ideal (x D - 22x D + 207x D - 975x*D + 1920, D , . . . 

05 : Ideal of R 

16 : time J = ideal Stafford I 

— used 73.08 seconds 

4 4 3 3 2 2 

06 = ideal (x D - 22x D + 207x D - 975x*D + 1920, 

3 15 2 15 2 14 14 13 

x D + x D + 15x D + 9x*D + 58x*D 

13 12 11 

+ 15D + 50D + D ) 

06 : Ideal of R 
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i7 : I == J 



o7 = true 



6 Conclusion 

The implementations of the algorithms constructed along the lines of the proofs 
of Theorems [| and ^ in Macaulay 2 work only on rather small examples for quite 
obvious reason: the expression swell in Grobner bases computations. 

Let us comment on the differences of algorithm of Hillebrand and Schmalc 
[|| and ours. Their algorithm takes care of (weaker) Theorem As a step it 
includes enumerating a certain infinite subset of polynomials in one variable and 
testing them to satisfy a certain property, where the testing procedure involves 
Grobner bases computations. Although we believe that their argument could 
be extended to build an algorithm for 100% of Stafford's theorem, it looks as 
the "test set" for the remaining 50% will be significantly more complicated. 
Hence, our constructive approach at every step of the algorithm seems to be 
more practical. Having programmed Hillebrand and Schmale's algorithm as 
well, we have to point out, that it faces the same type of expression swell as our 
program, hence the comparison of performance is just a theoretical question at 
this point. 

Finally, let us mention that the algorithm for finding a cyclic generator of a 
holonomic module is already included in the D-modules package for Macaulay 
2 H; eventually, the algorithms for finding two generators of a A„-ideal will be 
added to the package as well. 
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